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Dark matter concentrations and a search for cores in Milky 
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ABSTRACT 

We investigate the mass distributions within eight classical Milky Way dwarf 
spheroidal galaxies (MW dSphs) using an equilibrium Jeans analysis and we compare 
our results to the mass distributions predicted for subhalos in dissipationless ACDM 
simulations. In order to match the dark matter density concentrations predicted, the 
stars in these galaxies must have a fairly significant tangential velocity dispersion 
anisotropy (/? ~ —1.5). For the limiting case of an isotropic velocity dispersion (j3 = 0), 
the classical MW dSphs predominantly prefer to live in halos that are less concentrated 
than ACDM predictions. We also investigate whether the dSphs prefer to live in halos 
with constant density cores in the limit of isotropic velocity dispersion. Interestingly, 
even in this limit, not all of the dSphs prefer large constant-density cores: the Sculptor 
dSph prefers a cusp while Carina, Draco and Leo I prefer cores. The other four dSphs 
do not show a statistically significant preference for either cuspy or cored profiles. 
Finally, we re-examine the hypothesis that the density profiles of these eight MW 
dSphs can be quantified by a common dark matter halo. 

Key words: Galactic dynamics, dwarf galaxies, dark matter 



1 INTRODUCTION 



The currently favored cold dark matter cosmology (ACDM) 
has had much success i n reproducing the la rge scale struc- 
ture of the universe (see lKomatsu et alj|201ll . and references 
therein). However, on smaller scales there are some possi- 
ble discrepancies (e.g. Flores fc Primack 1994; Moorelll994l; 



2*1 



Klypin et alii 19991 ; iMoore et all 19991 ; iBovlan-KoIchin et all 
2011T ). One robust prediction from dissipationless ACDM 



simulations is that the inner density profile of dark matter 
halos shoul d be fairly cuspy at small radio with p(r) oc r~ 7 
and 7 ~ 1 dDubinski fc Carlberglll99ll;lNavarro et al.lll997l ; 
iDiemand et al.ll2008l ; ISpringel et alj2008f ). In contrast, some 
observed rotation curves from dark matter do minated galax- 
ies p r efer fits with 7 closer to than 1 (e.g. | de Blok et al.l 
200 ll; ISimon et alj|2003t 120051 ; iKuzio de Narav et al.ll200g[ 
de Blokll20101 ). Even when fits to the rotation curves are 
forced to 7 = 1 models, the resulting profiles are often less 
dense and less centrally-concentrated th an predicted from 
dissipantionless ACDM simulations (e.g.. lAlam et al.ll2002l ). 



At face value, the central density problem and cusp-core 
problem provides a challenge to cold dark matter theory. 

Milky Way dwarf spheroidal galaxies (MW dSphs) pro- 
vide particularly interesting dark matter laboratories. They 
include the most dark matter dominated systems known and 
are close enough to enable star-by-star kinematic studies 
(e.g. jWalker et alll2007l;ISimon fc G cha 200 jj; IStrigari et all 



200S 



Geha et al.l 120091; Ihokasi 12009): iKalirai et all 
tinez et all l201ll ; ISimon et all J201ll ; IWillman 



Willman et all l2010l : l Wolf et al.ll2010l . hereafter W10). Un- 



1 Interestingly, iPenarrubia et al.l 1I2OIO1 ) showed that dSphs are 
more likely to survive interactions with the MW disk if they are 
cusped rather than cored. 



fortunately, because these objects are dispersion-supported 
their unknown stellar velocity anisotropy makes it difficult 
to determine their mass distributions using a simple pro- 
jected Jeans analysis (e.g. IStrigari et alj|2007al . W10). 

In this paper, we ask whether a simple class of equilib- 
rium models prefer solutions with dark matter densities that 
are consistent with ACDM expectations, and we also explore 
whether the solutions allow constant dark matter density in 
the inner parts, in conflict with naive expectations of ACDM 
theory. 

Our work is mo tivated in part by the earlier study by 
iGilmore et al.l (|2007l ) (hereafter G07) where the results of an 
isotropic spherical Jeans analysis of the dispersion profiles of 
six of the eight brightest MW dSphs - Carina, Draco, Leo I, 



© 0000 RAS 



2 Wolf & Bullock 



Leo II, Sextans, and Ursa Minor - were presented. The main 
result of G07 was that solutions to the Jeans equation for 
these systems prefer dark matter models with an "approxi- 
mately constant density out to some break radius". In this 
paper, we include all eight bright MW dSphs (excluding the 
tidally-disrupting Sagittarius dSph). Our results are differ- 
ent and do not lend support to the claim that all MW dSphs 
inhabit halos with large constant density cores. We look at 
the full probability distribution of the inner slope of the 
dark matter density profiles and find the spread in allowed 
values to be large enough that we cannot make a definitive 
statement for half of these dSphs. Among the other half, our 
analysis shows that Sculptor strongly prefers a cuspy halo. 

We pay careful attention to the degeneracy of the dark 
matter density profile slope and the inner density profile of 
stars in our analysis. The inner density profile of the tracer 
stars (i.e. those with velocity measurements), as well as the 
stellar density profile derived from photometry, are uncer- 
tain in the inner parts. This uncertainty introduces signif- 
icant uncertainty in the determination of th e dark matter 
densit y profile - a point highlighted recently in lStrigari et al.l 
(|2010f ) (hereafter S10). We note that this degeneracy is sig- 
nificant, even in comparison to the well-known degeneracy 
of the mass profile with the unknown velocity dispersion 
anisotropy. 

Our conclusions are conservative in the sense that we 
force the inner light profiles to follow King (1962) models, 
which tend to bend sharply to cored light distributions as 
r — > and thus favor more core- like mass densities (see §2). 
King profiles are not unique fits for these systems, and by 
allowing a broader range of stellar profile shapes a wider 
range of density profile slopes become possible (S10). 

If there are constant density cores at the cen- 
ters of small dark matter halos it could signal the 
existence of non-st an dard da rk matter partic les (e.g . 
Hogan fc Dalcantonl |200d: lAvila-Reese et all |20Q]J ; 



Kaplinghatl 120051 ; Istrigari et all l2007d ). A 



modified 

dark matter solution would be particularly attractive if 
these cores were ubiquitous and of approximately the 
same density. Density cores might alternatively arise from 
baryonic physics that alters the dark matter distributions 
compared t o dissipationless simulations (as demonstrated 
explicitly by Mashchcnk o et al. 20081 ; iGovernato et al.ll2010] ; 



iPasetto et al.ll2010l ; lOh et alJ^Olll TTrie feedback processes 
that give rise to such cores rely on potentially stochastic 
details like the fraction of baryons that have been converted 
into stars (e.g. iBrook et al.ll2011f ). Therefore, if the cores 
are created by feedback processes, one might expect that 
the central density structure of dark matter halos will 
demonstrate significant variations from halo to halo, as 
seems to be the case for l ow surface brightness galaxies 
JKuzio de Narav et al1l2010l ). 

Lastly, we further explore the anisotropy-slope degen- 
eracy in order to address the question of what anisotropics 
are needed in order for simulations to reproduce these MW 
dSphs, if they are to live in ACDM-like halos. 

In the next section we provide a discussion of the de- 
generacies inherent in the use of Jeans modeling to constrain 
density profile slopes, even under the assumption of isotropy. 
We present our analysis and results in §3, where we discuss 
the mass modeling for each dSph, the question of a common 
halo for all the dSphs, the most probable stellar velocity 



dispersion anisotropies required to fit NFW profiles to the 
dSphs and a comparison of concentrations to those seen in 
recent ACDM simulations. We also discuss both the require- 
ments for our cuspy density profile fits to be physical when 
isotropy is assumed and the complexities in the stellar pop- 
ulation that we do not take into account in this work. We 
summarize our results in 54. 



2 EXPECTATIONS FROM THE JEANS 
EQUATION 

The spherical Jeans Equation relates the integrated mass of 
a spherically symmetric, dispersion-supported, collisionless 
system to its tracer velocity dispersion and tracer number 
density n+ir), under the assumption of dynamical equilib- 



M{r) = Lg. ( 7 , + 7<r -2£). 



(1) 



Here a r (r) is the radial velocity dispersion of the tracers 
and /3(r) = 1 — at /<*r quantifies the tangential velocity dis- 
persion. The additional terms quantify radial gradients in 
the stellar distribution 7* = — dlnn*/dlnr and velocity 
dispersion profile 7 CT = — dlna^/dlnr. While f3, and to a 
lesser extent y a , are difficult to determine empirically, 7* 
is in principle constrained by observations. If we assume a 
form of the anisotropy the n 7^ can be solved for directly 
(see iMamon fc Bouej|201ol , and Equations A7-9 in W10), 
and thus the Jeans equation immediately provides a path 
towards direct mass constraints. If we further assume that 
the stellar system is isotropic (/3 = 0), then our analysis 
simplifies considerably. 

However, even with the assumption of isotropy, con- 
straints on the slope of the mass density profile 7 = 
— d In p/A In r will require knowledge of the second derivates 
of both ni,(r) and oy(r), which will necessitate exquisite data 
sets. Typically, n*(r) is assumed to follow an empirically- 
motivated functional form, without significant concern for 
the allowed freedom of its detailed curvature at small ra- 
dius. For purposes of com parison , we w ill follow G07 and 
assume that n*(r) tracks a iKind (1 19621 ) profile. As we will 
now demonstrate, this is conservative in the sense that it 
forces a sharply flattening inner slope for the tracer stars, 
which in turn favors a flatter density slope 7 at small r than 
might otherwise be allowed. Indeed, S10 found that most of 
the available photometric data for the classical dSphs can 
be fit with divergent 3D number densities that asymptote 
slowly to power-laws. In this sense, the inner density slopes 
derived in this paper are biased to favor more core-like be- 
havior. 

One can understand how the shape of n* (r) affects core 
vs. cusp determinations by considering a proxy for the den- 
sity profile slope 7 ~ K = 3 — din M/dlnr. For /3 — 0, 
Equation [T] implies K = 2 + j a — (7* + 7^/(7* +7 CT ), where 
the primes indicate derivatives with respect to lnr. If we 
define 



i'*h* 



(2) 



in order to quantify the curvature of the light profile, we can 
rewrite the density slope proxy as 

k = 2 + la - 7 ; (7, + la y x - «* (1 + ja/^r 1 . (3) 
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Table 1 : Derived parameters of MW dSphs when modeled under the assumption of isotropy. 



Galaxy 



r 1/2 [pc] Cking «*(r 1/2 ) 7(r -+ pc) 7(100 pc) 7<r 1/2 ) 



5l0(Poore r oors M W) 



[PC] 



Carina 


334 


3.3 


1.9 


Draco 


291 


5.9 


1.2 


Fornax 


944 


5.2 


0.8 


Leo I 


388 


1.9 


2.7 


Leo II 


233 


3.5 


2.2 


Sculptor 


375 


13 


0.5 


Sextans 


1019 


9.6 


2.4 


Ursa Minor 


588 


1.1 


1.1 



00 +0.5 


0,1 


0.0+ 04 


0.3 


n 7 +0.3 


0.8 


0.0+ 04 


0.2 


5 +0 - 6 
u -°-0.3 


0.8 


1 1+ 01 
i - i -0.5 


1.2 


°-?±° .l 


0.8 


0.3+° n « 


0.5 



+0.3 
-0.2 
+0.3 
-0.2 
+0.2 
-0.4 
+0.3 
-0.1 
+0.5 
-0.4 
+0.1 
-0.4 
+0.4 
-0.3 
+0.4 
-0.3 



1+0.4 
-0.3 

HO. 3 

-ll. 1 
,1+0.3 

-0.2 
4+0.4 

-0.2 
,+1.6 

-0.5 
,+0.2 

-0.1 
i r+0.5 
1 -°-0.3 
-,+0.5 

-0.4 



0.4Z 

o.5~> 

1.4J 
0.4^ 
1^ 
1.3- 1 



i.o: 



1.2 



+ 1.2 



2.0 
1.1 
2.3 
1,1 
0.61 
0.61 
1.5 



o.i 

+0.4 

0.3 
+0.04 
-0.03 
+0.5 

0.6 
+ 1.7 

0.1 
+0.27 

0.04 
+0.27 

0.15 
+0.7 

0.1 



790+ 130 
,yu -110 



410 1 



-,+80 
-70 

DDU -200 



Note: Mean values of deprojected half-light radii (r 1/2 ) are from Table 1 of W10 while the mean King concentrations 

Cicing = ( r n m / r Kin E ) anc ' K * vam es are derived from the photometric properties listed in Table 1 of W10. Here K*(r 1 , 9 ) quantifies the 

curvature of the light profile at r 1 , 2 , as defined by Equation[2] As discussed in the text, galaxies whose light profiles have smaller values 

of the curvature will tend to favor cuspier overall density profiles. The total density log-slopes "f(r) = — dlnp/dlnr in the next three 

columns are derived when utilizing Equation [4] The density parameters in the last two columns are derived by utilizing the Burkcrt 

profile (Equation [5). Errors correspond to 68% likelihood values. 



At radii near the the deprojected half-light radius r ~ ri/2, 
the term multiplying k* in Equation [3] is positive and of or- 
der unitjQ which implies that a larger k* will favor a smaller 
N and a more core-like distribution. In the limit where 7^ 
and 7^. are small, k* will dictate the inferred density slope. 

To illustrate this point, consider the light profile of 
the Carina dSph, which prefers the parameter combina- 
tion r lim /r Kin = 3.27 with a deprojected half-light radius 
r i/2 = 334 pc when fit with a King profile (see Table 1 
of W10). S10 showed that the same photometry can be 
fit by a stellar light profile that is divergent at the center 
(with approximately the same half-light radius). At r ~ r 1/2 
both the S10 fit and the King fit have similar curvature: 
K*(r 1/2 ) = 1.05 and 1.12, respectively. But at the slightly 
smaller radus r = 0.5r 1/2 the light profile favored by S10 
has a significantly smaller curvature k+ — 1.5 compared to 
the King fit k* = 2.4, and this difference further increases 
at smaller radii, where k* diverges at small r for the King 
profile but drops to zero for the S10 fit. Indeed, compared to 
other common light distributions, the King fit imposes the 
strongest curvature towards small r and this will act to push 
solutions to wards mo r e core -like mass density profiles. For 
exam ple, a IPlummerl (119111 ) profile has K*(0.5r 1/9 ) = 1.4 
and a lHernquistl ( 19901 ) profile has K*(0.5r 1/2 ) = 0.3. Thus, 
for a fixed ov, King profiles will prefer smaller overall density 
profile slopes than all of the common alternatives. 

The above discussion highlights the difficulties asso- 
ciated with accurately constraining central density profile 
slopes using Jeans analysis. A similar point was emphasized 
by S10. In practice, one would need to know k* and j' a very 
accurately at small radii in order to definitively constrain the 
asymptotic density profile slope, even with the assumption 
of isotropy. Our approach instead is to adopt a conservative 
(and standard) assumption of the King tracer distribution, 
which will allow us to directly compare with the work of 
G07. We show that even in this limit, not all of the dSph 
galaxies prefer core-like distributions. 



This is because most commonly-used stellar tracer distributions 
yield 7* ~ 3 at r 1 / 2 (see Appendix B of W10) and empirically 
we expect |7 CT | < 1 owing to the rather flat velocity dispersion 
profiles observed (see Equations A7-9 of W10). 



3 DYNAMICAL ANALYSIS 

We use a spherical Jeans analyses to constrain the den- 
sity profile slopes of the eight classical MW dSphs: Carina, 
Draco, Fornax, Leo I, Leo II, Sculptor, Sextans, and Ursa 
Minor. The kinematic and photometric data we use are ref- 
erenced in Figure 2 and Table 1 of W10 jMunoz et al.ll2005l ; 
iKoch et al.l 120071 ; iMateo et all 120081 : IWalker et al.l l2009al ). 
We constrain the allowed mass profi le for each dwar f using 
a Bayesian analysis as described in iMartinez et al.l (|2009i ) 
and W10. 

We marginalize over a generalized density profile de- 
scribed by six parameters 



Ptot(r) = 



p s e 



(4) 



(r/r 3 ) c [l + (r/r\,) a ]( 6 - c )/ a ' 

using the same prior ranges discussed in Sec tion 2.3 o f W10 . 
When noted, we have also considered the iBurkertl (|1995f ) 
profile 



Ptot(r) 



Poo 



(l + r-Aco r o)(l + (r/r coro ) 2 ) 



(5) 



for comparison, with the same above priors for the scale radii 
and densities. In all of our analyses we marginalize over the 
photometric uncertainties of the fitted Kind (|1962j ) profiles 
(see Table 1 of W10). Finally, unless otherwise stated, we 
force the velocity dispersion tensor to be isotropic (j3 — 0). 

3.1 Density profiles for individual dSphs 

Table 1 summarizes results from our isotropic Jeans analy- 
sis for each of the eight galaxies we consider. The first three 
columns list relevant observational parameters, while the 
middle three columns list the derived constraint on the local 
density profile slope at three example radii: an asymptotic 
inner radius (r — > pc), an intermediate radius (100 pc), 
and at r . The last two columns display results from our 
analysis of constraining the density profile by Equation [5] 
As is shown in the table, and as is highlighted in our discus- 
sions below, some of the galaxies prefer cusps, some prefer 
cores, and others are indeterminate. Also, as expected from 
the discussion above, galaxies with the smallest curvature 
in their light profiles «;* (column 3) tend to be those that 
prefer the most cuspy overall density profiles (7 ~ 1). 

Posterior distributions for the asymptotic inner slope 
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Figure 1. Probability densities for 7(r) = — dlnp/dlnr as r — ¥ when modeling with Equation [4] the four MW dSphs with the highest 
quality kinematic data which is publicly available: Carina (black dotted), Fornax (red dashed), Sculptor (blue solid), and Sextans (orange 
dot-dashed), (a): Isotropic analysis (/3 = 0). This result shows that not all of the classical MW dSphs prefer to be cored (7 ~ 0) under 
the assumption that /} = 0. (b) and (c): Anisotropic analyses demonstrating how the likelihood for 7 widens as the range on the allowed 
anisotropics widens. 



are presented in Figure [T] for the four d Sphs that contain th e 
largest published kinematic datasets (| Walker et alj |2009a| ) ■ 
The left panel shows slope distributions for our isotropic 
analysis. Carina (black dotted) demonstrates a strong pref- 
erence for living in a cored dark matter halo (7 ~ 0), while 
Fornax (red dashed) and Sextans (orange dot-dashed) pre- 
fer moderate cusps (7 ~ 0.7), but their distributions are too 
wide to draw any strong conclusions. 

The most interesting case is Sculptor (solid blue) , which 
shows a strong preference for living in a cuspy dark mat- 
ter halo with 7 ~ 1. Sculptor's line-of-sight velocity disper- 
sion profile rises more rapidly in its center than any of the 
other well-studied MW dSphs, which drives the data to pre- 
fer cuspy profiles when isotropy is assumed. The physical 
reasoning for this is discussed in Section 3.1 of W10. Briefly, 
as one observes closer to the center of a system, the intrinsic 
radial dispersion projects onto the line-of-sight dispersion. 
Thus, a sharp rise in the observed dispersion could only be 
produced by two effects: the presence of a rising density (i.e. 
a cusp) or a highly radial stellar anisotropy. Since we are as- 
suming isotropy in our analysis, a cusp naturally arises. 

The middle and right panels of Figure [T] show how the 
asymptotic density slope distributions widen as we allow for 
wider ranges of (constant) velocity dispersion anisotropy, 
-1.0 < /3 < 0.5 and -3.0 < P < 0.75. Though not shown 
here, we also explored a case where the concentrated stellar 
component of Fornax was included in the total mass pro- 
file with a constant V-band mass-to-light ratio of T = 2. 
Despite the relatively low dynamical mass-to-light ratio of 
this galaxy (T^ ~ 9; see Table 1 of W10) we did not find 
any statistically significant change in the inner dark matter 
slope distribution. 

Figure[T]shows that under the assumption of isotropy at 
least one MW dSph (Sculptor) strongly prefers being hosted 
in a cuspy dark matter halo|j Inspection of Table 1 further 
reveals that four of the eight classical dSphs are consistent 
(one sigma) with inner cusps of 7 = 1, while three of the 
eight dSphs prefer asymptotic slopes of 7 ~ at one sigma. 

We also ran an alternate analysis where we set the mass 



Note that G07 did not analyze Fornax or Sculptor. 



density to be a Burkert profile (Equ ationQ along the lines of 
the analysis in lDonato et al.l (j2009f ) (hereafter D09) where it 
was proposed that over a wide range of luminosities, galaxies 
(including the six dSphs in G07) can be fit by cored Burk- 
ert profiles that have core densities and core radii related 
via p corc ex r~ r C such that p corc r corc = io 215±0 - 2 M pc -2 , 
independent of luminosity. Our results (see Table 1) show 
about a factor of four spread in the mean logio (p coro r coro ) 
values with the actual values for logio (p corc r corc ) about a 
factor of two smaller than those in D09. While the prod- 
uct of pcoro r corc is always well constrained, only three dSphs 
(Fornax, Sculptor, and Sextans) have constraining posterior 
distributions for r coro . 

Before moving on, we point out that each of these three 
galaxies have r coro ~ r 1/2 . Given that 7(r coro ) = 1.5 for a 
Burkert profile, these results do not indicate disagreement 
with dissipationless ACDM simulations, which show that 
dark matter halos have similar values for 7 around the same 
physical radii. 



3.2 Global population of MW dSphs 

We now move on to Figure [2] where the points presented 
in the left panel show the average density of each dSph 
within its deprojected half-light radius (p)i/2 as a function of 
r 1/2 . Our motivation for examining the density within r 1/2 is 
that constraints on the mass density within the deprojected 
half-light radius are less sensitive to assumptions about the 
velocity dispersion anisotropy, and thus they are more di- 
rectly constrained by data (W10). The right panel of Figure 
f2] shows the derived local log-slope at r 1/2 as a function of 
r 1/2 for the same set of galaxies under the assumption of 
isotropy. 

Recently, there have been several papers pointing out 
that the MW dSphs are consistent with inhabiting similar 
cuspy dark matter halos (jStrigari et al.ll2008l ; I Walker et al.1 
l2009bl . W10)Q Additionally, the mass within r has been 



4 Note that iMateo et al.l IU993T ) suggested that the classical 
dSphs had similar masses within their stellar extent. However, 
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Figure 2. Can a common halo reproduce all of the data? The solid black lines in each panel show the implied trends for a typical field 
NFW halo with maximum circular velocity V ma x = 26 km s _1 with r max = 5.7 kpc (the same NFW halo presented in Figure 3 of W10). 
We note that this is less dense than a typical subhalo of the same V max in recent dark-matter-only simulations of halos similar to the 
Milky Way. The two colored lines show cored Burkert profiles that obey the D09 constant surface density constraint with core radii 
of 170 pc (red dashed) and 850 pc (blue dot-dashed). Of these two, only the small cored halo provides a reasonable characterization 
of all the data in the left panel, but the same halo fails to reproduce the log-slope information in the right panel for any of the data. 
Meanwhile, the large cored halo only fares well at describing some of the data in the right panel. Of the three plotted, the cuspy NFW 
halo fares best at describing the full population within a common halo. However it is only consistent at the 2 — a level with predictions 
of the r — Vmax relation from recent simulations. 



shown to increase roughly as M 1/2 ~ rj (|Tollerud et al.l 
l201ll : IWalker et a"i]l20ld . W10), which is what one would 
expect if all of the galaxies were embedded within simi- 
lar dark matter halos that each obeyed the density scaling 



p(r) 



This is illustrated in the left panel of Figure 



[2] where the solid black line shows a NFW density pr ofile 
(which scales as r^ 1 at small r; see lNavarro et al.lll997h for 
a halo of virial mass 3 x 10 9 Mq and a typical concentration 
for a field halo (the same halo shown in Figure 3 of W10, 
which was chosen to represent all of the MW dSphs, not just 
the classical MW dSphs) . 

To more accurately assess the hypothesis that all of the 
MW dSphs lie in similar dark matter halos, one needs to 
simultaneously examine both the density and the log-slope 
of the population. We begin by noting the colored lines in 
Figure [2] which show the implied (p)(< r) vs. r and 7(7") vs. 
r relations for two example cored Burkert profiles - one with 
a small core (red dashed with r coro = 170 pc) and another 
with a large core (blue dot-dashed with r coro = 850 pc). The 
core densities for each of the Burkert profiles were chosen to 
follow the previously mentioned trend advocated by D09. 

As can be seen in the left panel of Figure [2] in order 
to match the density vs. radius relation for the sample of 
dSphs we need a Burkert profile that has a core radius that 
is smaller than the smallest half-light radius in our sam- 



the mass estimator used was not accurate (see Appendix C of 
W10 for an explanation), and thus the claim is not true. 



pie r core < 2 pc. A lthough this point was mentioned by 
IWalker et al.l (|2009bl ). if we were to include the ultrafaint 
MW dSphs as they did, the core size would need to be ap- 
proximately the size of Segue 1 (r 1/2 ~ 40 pc; see Table 1 of 
W10) in order to provide a decent fit to the data. Moreover, 
the same small-core Burkert that best fits the data in the 
left panel does poorly in the right panel because it is too 
steep to characterize the population. The large-core Burkert 
profile (blue dot-dashed) does better at fitting the data in 
the right panel but does much worse in the left panel because 
the implied core is too large. Overall, while it is possible to 
fit each individual dSph galaxy with a Burkert profile obey- 
ing a scaling of p C orc oc r<^, re (see Table 1), it is not possible 
to find a common cored halo solution that adequately re- 
produces the population (i.e. choosing Burkert profiles with 
intermediate core sizes fail in both plots). 

However, although it is certainly not a perfect fit to all 
of the data points in the right panel, the W10 NFW halo 
has a local log-slope profile that is consistent at one sigma 
with five of the eight dSphs while also faring well in the left 
panel (interestingly, six of the eight galaxies are one sigma 
consistent with a log-slope at r 1/2 steeper than 7 = 1). Also 
note that a NFW halo with a s cale radius smaller th an a kpc 
(like the one advocated for by Walker et al.ll2009bl ) will not 
provide a good fit to the data in the right panel of Figure [2] 

It should be pointed out, however, that the NFW profile 
plotted in this figure is lower than the median predictions 
from the Aquarius and Via Lactea 2 simulations. Thus, Fig- 
ure [2] makes it clear that /3 = is not a good assumption 
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when attempting to place MW dSphs in cuspy dark matter 
halos. We discuss this further in Section [3, 51 



3.3 Are the solutions physical? 

The solutions obtained by imposing the Jeans equation are 
not guaranteed to have a positive definite distribution func- 
tion. Here we investigate how a cuspy dark matter density 
profile is co nsistent with / 3 = and J*(r — » 0) — > 0. P revi- 
ous work by I An fc Evans! ||2006J ) [also see I White! (|l98lh and 
Ide Bruiine et al.l (|1996l ) for scale-free potential and stellar 
profiles] has shown that the log-slope of the external poten- 
tial should satisfy an inequality in order for th e distribution 
functi on to be positive definite. The results in lAn fc Evansl 
IJ2006T ) were noted in the limit r — > 0, but are applicable 
everywhere. Here we briefly recall their a rguments. 

We assume, as in I An fc Evansl (|2006r ). that the distribu- 
tion function of the stars has the form f(E, L) — L^ 2 ^ g{E) 
and that the dark matter provides an external potential for 
this stellar system such that the mass profile is monotoni- 
cally decreasing with radius. Given the assumption of a con- 
stant /3, 



n*(r)ocr- 2 ' 3 (-( ? (>(r)) 1/2 -' 9 / (l-E/<f>(r)) 



1/2-jS, 



(E)dE, 



(6) 

where <f>(r) is the gravitational potential sourced by the dark 
matter density profile, and the constant of proportionality is 
a positive definite function of /3. We neglect the influence of 
stars on the gravitational potential, which is a good approx- 
imation for the dSphs. With this, it is simple to show that 
k(r) — 7i*(r)r 2 ^(— ^(r))^ -1 ' 2 is a monotonically decreasing 
function of (f> for j3 < 1/2. Since d</>(r)/dr > and k(r) > 0, 
we have dlnfc(r)/dlnr < 0. This constraint implies: 



7*(r)> 2/3 + (1/2 -P)8(r), 



(7) 



where —S(r) — dln(— <j>(r))/A\nr > is the log-slope of 
the external (i.e. dark matter) gravitational potential. For 
j3 = 0, we have S(r) < 27* (r). For dark matter density 
profiles that diverge more slowly than 1/r 2 as r — > 0, we 
have S(r) — > 0. However, this is not sufficient since we have 
also assumed that 7*(r) — > as r — s- 0. 

The log-slope of the potential may be written as: 



5(r) 



V? (r) + AwG X°° p(x)x dx 
V c 2 (r) 



(8) 



4.5VJ 



0.5(r/r 



n2-3 







where V c (r) is the circular velocity and V max is the maximum 
circular velocity. To derive the final approximation on the 
second line we have assumed that the dark matter density 
profile diverges as (r/r a )~ for r -C r s . The constants in 
the approximations in the second line are a fit for a NFW 
profile but work well for all commonly used cored and cuspy 
dark matter halo density profiles. Since 7*(r) ex (r/rxing) 2 
as r <C r'King, we see that the constraint S(r) < 2j ic (r) will be 
violated at small radii because 27* for our tracer density (a 
King profile) asymptotes to more quickly than <5 as r — > 0. 
We may estimate the radius at which this constraint is 
violated by noting that 7*(r) ~ 3(r/rKm g ) 2 for typically ob- 
served values of the concentration of MW dSphs, and hence 
S(r) > 27* (r) for r < r s (0.3r-Kin g /r s ) 2 ' ■ For typical values, 



such as 1 — 1, rKing = 0.3 kpc, r s = 1 kpc, the constraint is 
violated at radii smaller than about 10 pc. 

Given that both photometric and kinematic data of 
MW dSphs are scarce in the region r < 10 pc (often the 
errors on the position of the center of the MW dSphs are 
larger than this) , substitution of a dark matter distribution 
with a small core or a profile that becomes more shal low with 
decre asing radius (such as the profile presented in Einastol 
Il965f ) would have little effect on our results outside of this 
region. In fact, recent high-resolution ACDM simulations do 
not lend support to the expectation that the inner density 
profile log-slope converges to -1, as an Einasto profile with a 
gradually decreasing log-slope amplitude seems to b e a bet- 
ter fit (e.g. lSpringel et aLllJOOJI ; iDiemand et al.ll2008l ). Thus, 
the implications of our isotropic results would not change 
were we to fulfill the requirements of a positive definite dis- 
tribution function. 



3.4 Multiple stellar populations 

Our analysis assumes that the entire spectroscopic sample 
follows a spatial distribution that is well described by the 
adopted photometric surface brightness fits. We have al- 
lowed for errors in the surface brightness fits, but not for 
the possibility that there are multiple populations with dis- 
tinct velocity dispersions. 

The motivation to model systems with multiple tracer 
components can be understood from the main analytic result 
of W10, which shows that the mass within a radius near the 
deprojected half-right radius can be determined nearly inde- 
pendent of stellar velocity dispersion anisotropy. If two sub- 
populations have different half-light radii, then one knows 
the mass well at two different locations, and thus the effects 
of the mass-anisotropy degeneracy are substantially reduced . 
This technique was first attempted bv lBattaglia et al.l (|2008l . 
hereafter B08) who modeled the Sculptor dSph assuming 
separate metal-rich and metal-poor stellar populations. In 
their analysis, B08 found that Sculptor prefers a cored mass 
distribution!^ 

The newer publicly available dataset for Sculptor 
(jWalker et alj|2007l . l2009al lch that we use contains approxi- 
mately three times more velocity measurements than in the 
B08 sample, and many new spectra are closer to the pro- 
jected center. We emphasize that this larger data set has 
not been used to study the issue of multiple populations. 
Further, B08 impose a strict prior of radial anisotropy (that 
increases with distance from the center) in their modeling 
and this could also skew the results in favor of cored profiles 
(see Section |3.5[) . 

The most important aspects to address when including 
information from multiple tracer populations in dynamical 
modeling include both deciding to which population each 
star (for which we have a velocity measurement) belongs and 
correctly modeling the underlying radial distributions of the 
subpopulations. For example, in B08 the two subpopulations 
are identified in the spectroscopic sample based on a sharp 



While this work w as being prepared for public ation we became 
aware of the paper bv lAmorisco &: Evans! a2012f) . which used the 
BOS data to perform a similar analysis. Our comments regarding 
BOS apply also to this work. 
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Figure 3. The posterior distribution of anisotropy (assumed to be spatially constant) for all of the MW dSphs considered in this paper 
when the density profile is fixed to be a NFW. The identification scheme is as follows - Carina: black dotted; Draco: red dashed; Fornax: 
blue solid; Leo I: orange dash dotted; Leo II: orange dotted; Sculptor: blue dashed; Sextans: red solid; and Ursa Minor: black dash dotted. 
The plot on the right shows that combined likelihood for /3 when each dSph is given equal weight. It is clear that the data demand a 
tangential anisotropy in the velocity dispersion tensor for cuspy profile fits. 



cut in metallicity (which tends to have a broad distribution 
in dSphs) and the radial density profile for each population 
is assumed to follow the radial distribution of the red and 
blue horizontal branch stars obtained through photometry. 
Uncertainties in this assignment translate to uncertainties in 
the estimated velocity dispersion profiles, especially that of 
the colder population. In addition, the radial distribution of 
the tracer populations should reflect these uncertainties (and 
they should not be completely fixed from fits to the photo- 
metric sample). Freedom in the radial distribution of the 
subpopulations would imply (as the analytic results of W10 
show) that the locations where the mass is well-constrained 
become more uncertain, thus making estimates of the mass 
profile and its derivative dependent on assumptions regard- 
ing the stellar velocity anisotropy. 

A related issue concerns foreground subtraction that 
may have benign effects for the entire population as a whole, 
but could affect one population more significantly than an- 
other. Further complications arise when one of the popula- 
tions is much more centrally concentrated, as is the case for 
the metal rich population identified in B08. The isophotal el- 
lipt icities could be diffe rent in the inner parts (for Sculptor, 
see lDemers et al.lll980T ) and adopting a half-light radius for 
each population without accounting for these effects could 
lead to a bias. A detailed investigation of these issues, fol- 
lowing the seminal work of B08, is required but it is beyond 
the scope of this paper. 



3.5 Stellar velocity dispersion anisotropy in a 
ACDM context 

Dissipationless ACDM simulations predict that dark matter 
halos generally follow cuspy density profile distributions over 
observable scales. For concreteness, we force the dynamical 



mass to be a NFW profile (a, — 1, b — 3, c = 1 in Equation 
31 and then derive the posterior distributions of f3 implied 
by the kinematic data. Our aim is to gain some insight on 
values of anisotropy required of these dSphs should they be 
embedded in a cuspy halo. 

Our results can be seen in Figure [3] where on the left- 
hand plot we superimpose the /3 distributions of all eight 
dSphs. The right-hand plot shows the combined likelihood 
of all weighted equally. For most of the dSphs, tangen- 
tially anisotropic velocity dispersions are favored. The 68% 
range for the combined sample is — 1 < /3 < 0. We note 
that this range is lower than what is generally reported for 
the anisotropy of dark matter particles, which typically are 
close to isotropic in the central regions and become radially 
anisotropic in the outer parts of the halo. However, there is 
no theoretical reason to expect the anisotropy of stars to be 
the same as that of dark matter. 

It is important to recognize that ACDM simulations 
predict not just the radial density profile of dark matter, 
but also the range of scale radii expected for a given subhalo 
mass or Vmax (i.e. the subhalo concentration). This amounts 
to a prediction for the log slope as a function of physical ra- 
dius. In Figure[3]we fixed the density profiles to be NFW but 
we did not force the scale radii to match those of subhalos 
predicted in ACDM simulations. Now we aim to investigate 
whether the derived 7(r 1/2 ) values are consistent with sub- 
halos of the appropriate concentration for their mass scale 
from high-resolution dissipationless ACDM simulations. 

The solid lines in each panel of Figure 2] show the pre- 
dicted log-slope run for a Vmax = 17 km/s Einsasto profile 
halo with a scale radius set by the median of the Aquarius 
simulation. Note that, even though it is not shown, this halo 
has the appropriate density normalization to provide a good 
match to the points in the left panel of Figure 2. It should be 
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Figure 4. Squares show our derived log-slope of the mass density at the deprojected half-light radius for each dSph assuming the general 
density profile in Equation [4] As marked, each panel presents a different choice of (tangentially biased) velocity dispersion anisotropy 
p. The solid black line, which is the same in each panel, shows the median expectation for V max = 17 km/s subhalos in the Aquarius 
simulation. The dashed lines show the one-sigma deviation from this relation. We see that in order for the MW dSphs to be consistent 
with the expectations for dissipationless ACDM subhalos, the stars must have a substantially larger tangential velocity dispersion than 
radial velocity dispersion. The purple circles, arbitrarily offset to the right for clarity, show a test that explores the sensitivity of our 
results to uncertainties in the central light profile. For these points, we have removed data for on-sky radii R < r 1 , 2 and we recomputed 
the log-slopes. We see that the overall result is insensitive to the inner data. 



noted that the radial dependence on the predicted log-slope 
from ACDM simulations does not change dramatically for 
subhalos that span reasonable V max ranges. 

The black open squares in the three panels of Figure 
2] show our derived results for j(t 1/2 ) for three choices of 
anisotropy (/? = —0.5, —1.0 and — 1.5). In deriving these 
points we have allowed a fully general halo profile (restoring 
freedom to all of the variables in Equation [J). As one moves 
from left to right in the panels, one can see that the only 
way for the dSphs to be consistent with the median subhalo 
expectation is for the velocity dispersion tensor of the stars 
to be tangentially anisotropic (/3 ~ —1.5). 

One concern with this result is that uncertainties in the 
stellar density profile shape will hinder our ability to accu- 
rately derive 7(r 1/2 ). The purple circles in Figure 4 (offset 
arbitrarily to the right for clarity) show results of a test de- 
signed to explore this possibility. Here we have removed data 
in the on-sky region R < |r 1/2 . We found that our results 
for 7(r 1/2 ) did not vary much (i.e. compare the purple cir- 
cles to the black squares, which were calculated with all of 
the data). This implies that the slope of the density profile 
at the deprojected half-light radius is determined predom- 
inantly by the local velocity measurements, and thus the 
uncertainty of the inner light profile does not contribute sig- 
nificantly to the derived values of 7(r 1/2 ) for the MW dSphs 
we are investigating. 

The main point to take from this subsection is that in 
order to decently model the classical MW dSphs as galaxies 
living embedded inside of subhalos in ACDM-based simula- 
tions (or even arbitrary NFW halos), a negative f3 for the 
stars is usually required. By imposing /3 = in the first half 
of this paper, we have been stacking the odds against finding 
cuspy halos. By simply adjusting all of the anisotropies to 
mildly negative values most of the dSphs can live in cuspy 
halos. Nevertheless, unless the negative anisotropy is signif- 
icant, those cuspy halos are typically less concentrated than 
expected for subhalos of the appropriate mass. In order to 



achieve an overall match to subhalo scaling relations, fairly 
significant tangential anisotropies are required (/? ~ —1.5). 



4 CONCLUSIONS 

We have performed dynamical Jeans analyses on the eight 
classical Milky Way dwarf spheroidal galaxies and we have 
explored their central density structure in comparison to 
predictions for subhalos in dissipationless ACDM simula- 
tions. We find that in the limiting case of an isotropic ve- 
locity dispersion (/3 — 0), the classical MW dSphs pre- 
dominantly prefer to live in halos that are less concen- 
trated than predicted in these simulations. In order to 
achieve an overall match to the cuspy profile shapes and 
subhalo concentrations expected from dissipationless ACDM 
simulations, fairly significant tangential velocity dispersion 
anisotropy is required (/3 ~ —1.5). This result is reminis- 
cent of the well-known cent ral density problem faced by 
ACDM (e.g. lAlam et alj|2002l ) in comparison to higher mass 
ga laxies. This effect is a ls o per h aps r elated to the findings 
of iBovlan-Kolchin et all (|201ll . I2012T ), who find that the 
classical MW dSphs are not dense enough to host the ex- 
pected number of massive subhalos around Milky Way size 
galaxies. The unexpectedly low concentrations of the dwarf 
spheroidal galaxies might be explained by stochastic feed- 
back effects, which in principle lower the densities of dark 
matter halos compar ed to what is predicted in dissipation- 
less simulations (e.g. [ Mashchenko et al.l 120081; iBrook et al 
20 111 ; iGovernato et al.ll2010l ; IPasetto et al.l l201dl : lOh et al 



20111 ). Another po ssibility is that the dark matter itself is 
non-standard fe.g. IVogelsberger et al.ll202ll ). 

Solutions to the isotropic Jeans equation show that not 
all of the MW dSphs prefer to live in halos that have con- 
stant density cores at small radii, even with a conservative 
treatment of the stellar density profile. As summarized in 
Table 1 and Figure [l] four of the eight dSphs are consistent 
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with 7 = 1 central cusps at one sigma. Sculptor strongly fa- 
vors a cuspy inner profile (7 ~ 1), while Fornax and Sextans 
prefer mild cusps (7 ~ 0.7), albeit with large uncertain- 
ties. Three galaxies (Carina, Draco, and Leo I) demonstrate 
a preference for a core (7 ~ 0). Therefore, our work lays 
to rest the possibility that the assumptions of isotropy and 
equilibrium force the inner density profiles of all of the MW 
dSphs to be cored. Rather, we have shown that the assump- 
tion of isotropy demands that their concentrations are lower 
than predicted in ACDM, regardless of the presence of a 
true core in the profile. 

Relaxing either the isotropy constraint or the shape of 
the inner stellar density profile broadens the allowed range 
of inner density slopes. If we adjust all of the anisotropics 
to negative values most of the dSphs prefer to live in cuspy 
halos. 

We also explored whether it was possible to find a 
common cored dark matter halo that simultaneously fits 
the average local density at r 1/9 and the local log-slope 
at r 1/2 under the assumption of isotropy. We showed that 
such a common halo is not possible when restricted to the 
family of Burkert profiles advocated by D09, but that the 
NFW profile used by W10 is a more reasonable fit. If we 
drop the log-slope constraint derived assuming isotropy, 
then we can find a single Burkert profile which matches the 
density-radius relation for the classical MW dSphs, as long 
as the core radius is smaller than the typical size of any 
galaxy in the population (r corc < 200 pc). Interestingly, one 
of the common Burkert halos with which we fit the local 
density at r 1/2 is actually steeper than the W10 NFW halo 
over the radii probed by the data. 
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